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Scanning tunneling microscopy and spectroscopy of sodium-chloride overlayers on the 
stepped Cu(311) surface: Experimental and theoretical study 
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The physical properties of ultrathin NaCl overlayers on the stepped Cu(311) surface have been 
characterized using scanning tunneling microscopy (STM) and spectroscopy, and density functional 
calculations. Simulations of STM images and differential conductance spectrum were based on the 
Tersoff-Hamann approximation for tunneling with corrections for the modified tunneling barrier 
at larger voltages and calculated Kohn-Sham states. Characteristic features observed in the STM 
images can be directly related to calculated electronic and geometric properties of the overlayers. The 
measured apparent barrier heights for the mono-, bi-, and trilayers of NaCl and the corresponding 
adsorption- induced changes in the work function, as obtained from the distance dependence of 
the tunneling current, are well reproduced by and understood from the calculated results. The 
measurements revealed a large reduction of the tunneling conductance in a wide voltage region, 
resembling a band gap. However, the simulated spectrum showed that only the onset at positive 
sample voltages may be viewed as a valence band edge, whereas the onset at negative voltages is 
caused by the drastic effect of the electric field from the tip on the tunneling barrier. 

PACS numbers: 68.35.Ct, 68.55.Ac, 68.37.Ef 
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I. INTRODUCTION 

Thin films of insulating materials on metal surfaces 
are of direct technological interest in microelectronics, 
catalysis and as a corrosive protection. They are also at- 
tracting increasing interest for nanoscience and nanotech- 
nology as potential substrates for atomic and molecular 
manipulation by the scanning tunneling microscope. In 
particular, this interest stems from the possibility to de- 
crease the electronic interaction of assembled molecules 
and adatoms with the metal support. The non- vanishing 
electron density that reaches through the ultrathin insu- 
lating films, thickness below 1 nm, allows for these sur- 
faces to be studied by means of STMi^S*^. A most im- 
portant prerequisite for the realization of this possibility 
is to be able to grow stable and atomically thin insu- 
lating films on metal surfaces with a well-characterized 
and ordered geometric structure. So far the number of 
systems studied that fulfills these conditions have been 
very limited^. One interesting candidate is provided by 
NaCl overlayers on a stepped Cu(311) surface, which we 
have studied in detail by scanning tunneling microscopy 
(STM) — and also most recently by density functional 
calculationsfi. 

These experimental and theoretical studies revealed 
several interesting aspects of the growth, structure and 
bonding of these systems. The initial growth of NaCl on 
the Cu(311) surface was found to be two dimensional and 
commensurate. The lattice match of the overlayer with 
the substrate was suggested in the experimental study to 
be stabilized by the incomplete screening of the step Cu 
atoms resulting in an electrostatic interaction between 
these atoms and the CI atoms. Density functional cal- 



culations corroborated this suggestion and showed also 
that the bonding of the overlayer on the surface was fur- 
ther stabilized by the formation of a weak chemical bond 
between the step Cu atoms and the CI stoms. These cal- 
culations revealed also large relaxations of the monolay- 
ers such as a buckling of the layer with a large influence 
on the work function. These results call for a detailed 
scanning tunneling microscopy study of the topography 
of and the tunneling conductance through these overlay- 
ers. In particular, the imaging of insulating overlayers 
on metal surfaces for biases well below 1 V has not been 
studied before. 

In this paper, we present atomically resolved topo- 
graphical images of NaCl mono-, bi- and tri-layers on 
Cu(311) obtained by an STM and discuss them in rela- 
tion to simulated images based on density functional cal- 
culations and the Tersoff-Hamann (TH) approximation^. 
In this approximation, STM images and scanning tunnel- 
ing spectra are obtained from the local density of states 
(LDOS) at the position of the tip apex. The conductiv- 
ity of the sample are further investigated using scanning 
tunneling spectroscopy. The measured 4^ spectra are 
compared with simulated spectra and discussed in rela- 
tion to the energy dependence of the calculated LDOS. 
This study includes also a calculation and discussion of 
the apparent heights and work functions of these layers 
obtained from tunneling measurements. The apparent 
height of an NaCl overlayer contains information about 
the "local" work function, that is, the tunneling barrier. 
The tunneling barriers are obtained from the measured 
distance dependence of the tunneling conductance. 

STM investigations of insulating overlayers provide im- 
portant information about the growth mechanisms and 
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FIG. 1: Experimental STM images and calculated LDOS im- 
ages of a NaCl monolayer in the p(3x2)-II structure, (a) Cal- 
culated LDOS image at an average distance of 7.9 A from the 
Cu(311) surface. The surface unit cell used in the calculation, 
see Fig. HI (right), is indicated, (b) Experimental STM image 
with tunneling parameters 210pA and -24meV. In both (a) 
and (b) the presented area is ISxlsA'^. 



FIG. 2: Experimental STM images and calculated LDOS im- 
ages of a NaCl bilayer on Cu(311). (a) Calculated LDOS 
image at an average distance of 9.0A from the Cu(311) sur- 
face. The surface unit cell used in the calculation, see Fig. |5] 
(left), is indicated, (b) Experimental STM image with tun- 
neling parameters 360pA and -lOOmeV. In both (a) and (b) 
the presented area is ISxlsA^. 



the local electronic structure of the over layer. In earlier 
STM investigations of ultrathin insulating overlayers on 
metal surfaces such as MgO on Ag(lOO)^ and NaCl on 
Al(lll) and Al(100)i the experiments were carried out 
at relatively high biases of about leV to provide atomic 
resolution of the overlayer. Atomically resolved STM im- 
ages have also been observed for NaCl layers on Ge(lOO)*, 
where only one type of atoms was found to be imaged 
as a protrusion"'^ . A combined STM and density func- 
tional study of NaCl layers on an Al(lll) and Al(lOO) 
surface revealed that only the CI atoms were imaged as 
protrusions. In this work, the atomic contrast as well as 
the apparent height of the NaCl layers was discussed by 
comparing the experimental images with calculations of 
LDOS within the TH approximation. MgO layers on an 
Ag(lOO) surface were studied using measured ^ spec- 
troscopy and the band gap of the bulk MgO electronic 
structure was found to be developed within the first three 
monolayers^. This finding was corroborated by density 
functional calculatrions of the energy dependence of the 
LDOS within the MgO layers. 

The paper is organized as follows. In section llll our 
experimental STM and ^ results for mono-, bi- and tri- 
layers of NaCl on Cu(311) are presented. In section llTTl we 
present the geometric and electronic structure and simu- 
lated STM images and M^ spectra for the same systems. 
The experimental and theoretical results are compared 
and discussed in Section Hvl Finally, section Ivl conchides 
the paper. 



II. EXPERIMENTAL METHODS AND 
RESULTS 

Our experiments were carried out with a low temper- 
ature scanning tunneling microscope^, operated at 13K. 
We used a chemomechanically polished Cu(311) single 
crystal, which was cleaned by Ne+ sputtering and an- 
nealing at 750K. NaCl was evaporated thermally and 
the deposition and annealing temperature was varied in 
the range of 400K-570K. In this temperature range Cu 
surface atoms are mobile. As an STM tip we used an 



electrochemically etched tungsten wire. Bias voltages re- 
fer to the sample voltage with respect to the tip. The 
STM images were taken at 13K. The STM images for 
mono- and bi-layers that are presented in Figs, ^b)- 
|2Ib) and bare Cu(311) were acquired at relatively low 
voltages of I V \< lOOmV, being well within the band 
gap of bulk NaCl. For the bare Cu(311) and the mono- 
layer structures our experimental images were obtained 
under identical tunneling conditions, that is, both ap- 
plied current and voltage, and tip were identical. Since 
the growth of the first layer was always completed before 
the formation of the second layer we could not be sure 
that the tip structure was identical in the imaging of the 
bare and the bilayer covered surfaces as for the bare and 
monolayer-covered surfaces. However, we could simulta- 
neously record STM images of islands with one to three 
layers on top of the monolayer. Up to four layers of NaCl 
could be imaged with atomic resolution. 

In our experiments we observe two different structures 
for the NaCl monolayer, shown in Figs. H^b) and|3|[b). 
As argued in our earlier study, we attribute these images 
to be associated with the structures presented in Fig. 
Here we follow the notation of Ref.^ and designate these 
structures as p(3x2)-II and p(3x2)-I, respectively^"'. In 
the STM images of these structures, we identify one char- 
acteristic feature for each structure. In the STM image 
of the p(3x2)-II structure every second protrusion appear 
slightly brighter. For the p(3x2)-I structure we note a 
slight tendency for the protrusions to group into pairs, 
even though, this effect is small and is just at the limit 
of the present resolution. In contrast to the monolayer, 
atomically resolved STM images of the bi- and trilay- 
ers (Fig. |2Ib)) did not show any differences between the 
domains corresponding to p(3x2)-I and -II structures, al- 
though it could be concluded that both these domains 
are present with the help of characteristic defects, con- 
sisting of a missing Cu atoms in the topmost substrate 
layer^. 

The apparent heights, Azi, of the NaCl overlayers 
(I — 1,2,3 for the mon-, bi-, and tri-layer) are presented 
in Table QJ Azi is defined as the vertical change in tip 
position when a step edge is scanned. We find that Az/ 




FIG. 3: Experimental STM images and calculated LDOS im- 
ages of a NaCl monolayer in the p(3x2)-I structure, (a) Cal- 
culated LDOS image at an average distance from the Cu(311) 
surface of 7.9 A. The surface unit cell used in the calculation, 
see Fig. |S](left), is indicated, (b) Experimental STM image 
with tunneling parameters 210pA and -24meV. In both (a) 
and (b) the presented area is 15xl5A^. 
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TABLE I: Calculated work function ^theory and experimen- 
tally measured mean barrier height $oxp. of NaCl layers on 
Cu(311). The values in brackets indicate the relative changes 
with respect to the clean copper surface. Azi are the apparent 
heights obtained from LDOS calculations (theory) and STM 
experiments (exp) and Azf""^' are the calculated geometric 
height. /=1,2 and 3 for mono-, bi- and trilayers, respectively. 
The result for the NaCl monolayer is for the p(3x2)-II struc- 
ture. 

Cu(311) mono bi tri 

fthcory (eV) 4:31 3.53 (-18%) 2.85 (-34%) 2.83 (-34%) 

'l>cxp.(eV) 3.81 3.23 (-15%) 2.82 (-26%) 2.82 (-26%) 



Azi (exp) (A) 
Azi (theory) (A) 



^^gcom. 



(A) 



1.6 
1.8 
2.5 



1.4 
1.1 
2.8 



1.3 



2.8 



decreases with / and is positive for all l. As also shown 
in Tabled the geometric heights, Azp°"\ of the NaCl 
overlayers are about lA larger than Azi. This implies 
that the tip-surface distance decrease with I. 

To quantify work function changes upon adsorption of 
NaCl layers, we have carried out measurements of the 
distance dependence of the tunneling current /^". The 
feed-back loop was switched off and a voltage ramp was 
applied to the z-piezo drive while measuring the current 
/ as a function of tip-sample distance z. The resulting 
sweep in z measured 2.7 A in height and had a duration of 
about one second. The resulting current curve followed 
closely an exponential dependence, I{z) oc exp(— 2kz) 
from which an apparent barrier height 4> was defined 
from the decay constant k as <& = Ti^k^ /{2me)- The 
relationship between $ and the work function $ of the 
sample is not straightforward and has been discussed at 
length in the literatur o^^i^^ . The systematic deviation 
of $ from $ decreases with increasing tip-surface dis- 
tance, corresponding to a larger vacuum region between 
the tip and the sample. So this deviation was minimized 
by recording spectra at a lowest possible current, rang- 
ing down to lOpA, so as to maximize the tip-surface dis- 
tance. Although $ can not be determined directly, the 
adsorption-induced relative changes in $ are expected 
to correspond to the relative changes in $. The exper- 



FIG. 4: Experimental dl/dV spectrum for a monolayer NaCl 
on Cu(311). The tip-sample distance corresponds to a current 
setpoint of / = 0.3 nA &t V = 3.8 V. This distance is by 
about 6 A larger than the distance for the imaging parameter 
set used in the Figs. 0and|21 



imental uncertainty in determining $ from k is set by 
the vertical length scale calibration, being a few percent. 
Using this procedure, we find that the measured work 
function difference between a NaCl monolayer and the 
clean surface agrees well with the work function differ- 
. ence measured by an atomic force microscopy tipi^. 

The apparent barrier height $ for the bare surface and 
NaCl overlayers deduced from the measured the I{z) are 
shown in Table HJ These results show that there is a 
substantial decrease in $ upon adsorption of the first 
layer and an even larger decrease for the second. No 
further changes of $ are observed for the thicker layers. 
This result indicates that the nature of the monolayer 
differ from the interface layer in the multilayers. 

The insulating properties of the monolayer have been 
investigated using dl/dV spectroscopy, shown in Fig. ^ 
for a monolayer. When acquiring the spectra the tip was 
retracted about 6A. The spectra is not sensitive to the tip 
position relative to the surface. The poor conductivity 
of the monolayer is manifested by a very small tunnel- 
ing current between -4 and -I-3V and exponential off and 
onsets of ^ at these voltages. The other characteristic 
feature in the spectra is the shoulder at about -t-3.2V. 



III. THEORETICAL METHODS AND RESULTS 

The results obtained from the STM measurements 
for the bare surface and the NaCl overlayers such as 
constant current images, work functions and apparent 
heights have been analyzed by extending our recent den- 
sity functional calculations of the electronic and geomet- 
ric structures of these system^, including now also the 
trilayer, to tunneling. These calculations were based on 
a plane wave basis set and the projector augmented wave 
(PAW) methodi^ii^ as implemented in the VASP^'^ code. 
A generalized gradient approximation^^ was employed for 



TABLE II: Calculated interface energies and geometrical pa- 
rameters of a NaCl monolayers, bilayer and trilayer on the 
(311) surface of Cu. The interface layer for the bilayer and 
trilayer correspond the p(3x2)-I structure. The adsorption 
energy, Ea, is given in eV per NaCl pair of the interface and 

is defined as Ea = {Eu^aCl/Cu — EnNaCl — E(i^n)NaCl/Cu)/T^, 

where 1 is the number of NaCl layers in the overlayer, m is 
the number of NaCl in the interface and n=l, n=l,2 and 
n=l,2,3, for monolayers, bilayers and trilayers, respectively. 
A& = zci — ZNa, where zci and zjva are the coordinates of 
Na and CI atoms perpendicular to the surface, zci is the dis- 
tance between the CI atom and the Cu surface layer, dci-ci 
is a Cl-Cl distance along a row in the layer. The distance 
between CI atoms on two neighbouring Cu-rows is 4.26 A. 
The superscripts (i) and (6) for the p(3x2)-I structure corre- 
sponds to CI atoms in top and bridge sites relative to the top 
Cu(311) layer, respectively. The superscript (i') and (b') for 
the p(3x2)-II structures corresponds to Na atoms in top and 
bridge sites relative to the second Cu(311) layer, respectively. 
See Fig. 1 in Refi^. 

Eg (eV) A b (A) ZCI (A) dci-ci 

monolayer 



p(3x2)-I 
p(3x2)-II 



0.42 
0.41 



0.31'' >, 0.35'*'' 2.54 3.78,3.94 

0.36<*\ 0.29'''' 2.56'*' 2.49*''' 3.86, 3.86 



bilayer 

interface 035 0.13(''>, O.::'*"'' 2.50('' 3.83, 3.89 
top 0.22 < 0.02 5.29 3.86, 3.86 



trilayer 

interface 034 0.15^''', 0.13^'''' 2^51 3.83, 3.89 

middle 0.25 -0.09 5.25 3.86, 3.86 

top 0.31 0.08 8.14 3.86, 3.86 



the exchange-correlation potential. The system was rep- 
resented by a slab in a super cell geometry^^ and the 
equihbrium geometry was obtained by a structural opti- 
mizatioEk2&. We considered the two p(3x2)-I and p(3x2)- 
II structures for the NaCl monolayer, shown in Fig. 
and the structures for the bilayer and trilayer that were 
formed by adding subsequent NaCl layers on the p(3x2)- 
I structure of the monolayer. In Tables Q] and |nl we 
show the calculated work functions and key bonding pa- 
rameters for the NaCl layers. For easy reference these 
tables include also the parameters for monolayer and the 
bilayer, previously presented in Ref.;^. The calculated 
partial s, p, and d density of states of the CI atoms in 



trilayer are shown in Fig. [7| As expected, the width of 
the CI peak is largest for the interface layer (Fig. [7|[c)) 
and smallest for the outermost layer (Fig. CJa)), due to 
larger overlap with the Cu states for the interface layer. 

The STM images and the apparent heights of the 
overlayers were simulated using the TersofF and Hamann 
(TH) approximation^ for the tunneling current. In this 
approximation the tunneling current for low bias is pro- 
portional to the LDOS, p{rQ,ep), of the surface at the 
position of the tip apex, Tq. The LDOS was calculated 
from the Kohn-Sham wavefunctions, see Ref>i^ for de- 
tails. The topographical STM images were obtained from 
the topography of the contours of constant /9(ro,eF), 
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FIG. 5: Ordered structures of a NaCl monolayer on Cu(311): 
(left) p(3x2)-I and (right) p(3x2)-II. The dark and light filled 
circles represent Cu atoms of the first and second layer, re- 
spectively, and the large white and small black open circles 
represent CI and Na atoms, respectively. The unit cells used 
in the calculations are indicated. In the p(3x2)-I structure all 
CI atom sites are equivalent, whereas there are two inequiv- 
alent Na atom sites, top (f) and bridge (b'), with respect to 
the second Cu layer. In the p(3x2)-II all Na atom sites are 
equivalent, whereas there are two inequivalent CI atom sites, 
top (t) and bridge (b) sites with respect to the first Cu layer, 
see also Fig. 1 in Ref,^. 



which henceforth will be referred to as LDOS images. 

We have also simulated differential conductance, ^ 
spectra, using the LDOS outside the surface. For these 
simulations it is necessary to consider the effects on the 
tunneling caused by the applied electric field from the tip. 
Following Langl^, these effects are partly accounted for 
by calculating -rrr within the average tunneling barrier 
approximation as 



dl d r^F+eV 

dV - dV^l Psieme,V)de) 

= epsiep + eV)T{eF + eV,V) + T -^r{Ps{e)T(e, V))de , 



(1) 
(2) 



r 



where T(e, V) is the transmission coefficient through the by 
average barrier for an electron with energy e and is given 



r(e, V) = exp {-2s^2m{^av + e^ - e + eV/2)/n) (3) 



CO 



CO 



o 

01 



> 



°6 



I I 
I I 
If 
I I 

V\ ' / 

: \ il 

• \ ll 

- \ '/ 

- \ 1/ 



-4 



-2 

Voltage(V) 



FIG. 6: Calculated -jv spectra for a monolayer NaCl in the 
p(3x2)-I structure on Cu(311). The solid and dotted lines are 
calculated from Eq. Q using the full energy dependence of 
ps(e) and a constant ps(e), ps(e) ~ Psi^p), respectively. The 
dashed-dotted line is the result obtained using only the local 
density of states contribution to ^ 
the calculations we used Zs = 2 A, 
eV. 



"^^ (first term in Eq. @). In 
s = 10 A and <E><i„ = 3.5 



Here <J>au = ($t + $s)/2 where <I>t and $5 is the tip and 
sample work function, respectively, ps(e) is the LDOS at 
a distance Zg just outside the surface layer defining the 
boundary of the vacuum region and s is the distance be- 
tween Zs and the tip apex. The first term in Eq. |2Jl is the 
LDOS at the tip-apex taking into account electric field 
in the tunneling junction. In the limit of vanishing volt- 
age only this term survives and the TH approximation is 
recovered. Note that this result for ^^ is well defined for 

dV 

\V\ < 2$ai, but the average barrier approximation breaks 
down for voltages approaching the field-emission regime, 
that is, for V > $s/e or V^ < -$T/e. 

In Fig. El we show the simulated ^ spectra for a 
monolayer of NaCl in the p(3x2)-II structure. The be- 
havior of ^ in the voltage region -3 to 3 V is well under- 
stood from the behaviour of the LDOS. In particular, the 
rapid increase of ^ around y = 2 V is an LDOS effect 
and may be viewed as conduction band edge. In con- 
trast, the rapid decrease of ^ around V — -4.5 V is not 
an LDOS effect associated with any valence band edge. 
Note that the large contribution of electron states of p or 
s character around the CI atoms to the partial DOS in 
Fig. [3 in the corresponding energy region does not show 
up in the LDOS. As illustrated by the result for ^ in 
Fig. El using an energy-independent ps(e) = psi^p), the 
rapid decrease comes from the second term in Eq. (|21l, 
which describes the drastic effect of the voltage depen- 
dent barrier on the tunneling and the larger probability 
for tunneling of electrons with energies closer to the vac- 
uum level at larger negative voltages. In the calculation 
we have used the same work function for the tip and the 
sample so we expect even a more dramatic onset of M^ for 
V < -3.5 V corresponding to the field-emission regime. 




FIG. 7: Partial s, p and d density of states around CI atoms 
in a NaCl trilayer on Cu(311). The CI atom is in (a) the third 
(outermost) NaCl layer, (b) the second layer, and (c) the first 
(interface) layer. Note that the scale of the presented PDOS 
on the three CI atoms differ by several orders of magnitude. 



IV. DISCUSSION 

We begin by comparing and discussing the measured 
apparent barrier heights and calculated work functions 
for the p(3x2)-I monolayer, the bilayer and the trilayer. 
As shown in Table JH the calculated changes in work func- 
tions reproduce nicely the measured changes for the bare 
surface and its change upon adsorption of the overlayers. 
Note that the measured apparent barrier height and the 
calculated work function can not be compared directly. 
From the calculations we are able to understand the re- 
duction of <& upon adsorption and its layer dependence in 
terms of bonding and buckling of the monolayer. This ex- 
planation was discussed in detail in our previous studji^ 
and will only be summarized here. 

In general, the surface contribution to <& is determined 
by the surface dipole layer, which have contributions both 
from electronic and ion-core rearrangements upon ad- 
sorption. The reduction of $ upon adsorption was ar- 
gued to be caused by the formation of a weak chemical 
bond between the CI atoms and the metal substrate and 
the screening of the Na+ ion by the metal electrons. The 
substantially smaller reduction of <I> for the monolayer 
than for the bilayer and the trilayer is caused by the large 
buckling of the monolayer with a large inward relaxation 
of the Na"*" ions relative to the Cl~ ions, resulting in a 
dipole layer that increases the work function. This buck- 
ling is substantially reduced for the bilayer and the tri- 
layer and the work function reductions for these layers are 
governed mostly by the electronic rearrangements, which 



are limited to the interface region. The measured <I> is 
thus consistent with a buckhng of the monolayer and an 
interface bonding that involves electron rearrangements 
between the overlayer and the substrate. 

Next we turn to a comparison and a discussion of 
observed and simulated STM images for the bare sur- 
face and the overlayers. These images contain informa- 
tion about the surface electronic and geometric structure 
through the tails of the wave functions in the vacuum 
region. Since the TH approximation for the tunneling is 
not based on any detailed model for the tip s density of 
states, we are not able to determine the distance, zq, be- 
tween the tip apex and the surface Cu layer from the tun- 
neling conductance. However, the bare surface and the 
two monolayers were imaged under identical tunneling 
conditions so we can use the image for the bare surface 
to calibrate the value of the /o(ro, e_F) in the calculation of 
the LDOS images and to make a prediction of the images 
for the monolayers. 

For the bare Cu surface, the calculated LDOS image 
gives the same tip-surface corrugation across the Cu rows 
as the STM image at an average zq of about 6.1 A. No 
corrugation could be resolved along the Cu rows, neither 
in the LDOS image nor in the STM image. In contrast to 
the flat single crystal metal surfaces for which the TH ap- 
proximation fails to reproduce atomically resolved STM 
imagesSS, the measured tip-surface corrugation can be 
reproduced by a realistic value for zg, for which the TH 
approximation should be applicable. For the flat metal 
surfaces, the imaging mechanism has been argued to in- 
volve a direct tip-surface interaction^i. 

In Figs. Hand 13 we show the calculated LDOS images 
for the investigated monolayers, using the same value of 
p(ro, ep) as for the bare surface that reproduced the ob- 
served surface corrugation. The corresponding average 
value of zq is 7.9 A. The resulting LDOS images are in 
overall good agreement with the observed STM images 
although the corrugation is somewhat larger in the LDOS 
images than in the STM images. The LDOS images sup- 
port our original interpretation of the observed STM im- 
ages that the protrusions correspond to CI atoms^. The 
dominant contribution from the states of the CI atom to 
the tunneling is understood from the participation of the 
p states in the bond as indicated by the corresponding 
partial density of states in Fig. [T] 

As shown in Fig. Q] the LDOS image for the p(3x2)-II 
monolayer reproduces nicely the experimental observa- 
tion that every second CI is imaged slightly brighter. In 
our calculations we attribute the brighter protrusions to 
be derived from CI atoms at bridge sites. This difference 
in brightness corresponds to a difference in zq of about 
0.054 A and is consistent with larger adsorption height, 
zci, of 0.07 A (see Table llT)l for the CI atoms at the bridge 
sites than at the top sites. However, the hybridization of 
the 3p states of the CI atoms with the metal states at the 
Fermi level is found to be stronger for the top site than 
for the bridge site. This electronic effect counteracts in 
part the geometric effect on the tip-surface corrugation. 



As shown in Fig. O the tendency of pairing of the pro- 
trusions in the STM image for the p(3x2)-I monolayer is 
well reproduced by the LDOS image. This pair formation 
reflects the geometric relaxations of the CI atoms along 
the Cu rows, resulting in two different Cl-Cl distances, 
which differ by 0.16 A (see Table HH . However, the ge- 
ometric relaxations contribute only a minor part to the 
difference in distance of about 0.9 A between protrusions 
along the rows. So in this case the dominant part of 
this difference is an electronic structure effect associated 
with the tilting of the Cu-Cl bond and the hybridized 3p 
states of CI with respect to the surface normal. 

In comparing STM and LDOS images of the bilayer 
we cannot adopt the same unambiguous procedure for 
the choice of average zq as we used for the monolayers. 
The bare surface and the bilayer are never imaged under 
the same tunneling conditions because the monolayer is 
fully developed before the formation of the bilayer. How- 
ever, in our comparison with the STM image we use the 
same value for p{ro,eF) as above, which corresponds to 
an average zq of 9.0 A. As for the monolayers we find that 
the protrusions in the LDOS image for the bilayer shown 
in Fig. 13 correspond to CI atoms. As shown in Fig. [7| 
there is still a contribution from CI states to the partial 
DOS around the Fermi level. This contribution is much 
smaller than for the CI atom in the interface layer but 
is compensated by the average distance of the tip from 
the outermost NaCl layer being only about 3.5 A for the 
bilayer, which is about 2 A shorter than the correspond- 
ing distance for the monolayer with the same value of 
p(ro,eF)- The agreement between the LDOS and STM 
image with respect to amplitude of the surface corruga- 
tion is about the same as for the two monolayers but in 
contrast to the STM image the protrusions are more dif- 
fuse and paired into dimers in the LDOS image. This 
pairing is a pure electronic structure effect because the 
distances between CI atoms across the rows are the same 
even though the atoms are inequivalent. Note that the 
distance of about 3.5 A between the tip and the outer- 
most NaCl layer is slightly too short to justify the TH 
approximation and an enhancement of the contrast from 
tip-surface interaction cannot be ruled out^^. In particu- 
lar, the tip-surface distance is even shorter for the tri- and 
quadrolayers so the observed atomic resaolution might be 
caused by the tip-surface interactions being sensitive to 
the identity of the underlying atom. 

The apparent heights Azi of the single steps on the 
NaCl multilayer I provide indirect information about the 
conductance through the layers. For example, Az/ is 
equal to zero or the geometrical height for a multilayer 
/ of vacuum or Cu, respectively. Because bulk NaCl is 
a wide-gap insulator with a calculated band gap of 5.0 
e\i21, the electrons around the Fermi level have to tun- 
nel through the film so its conductance should be sub- 
stantially less than an Cu layer, resulting in a smaller 
apparent height than the geometrical height. 

As shown in Table Q] the calculated Az/ are in good 
agreement with the experimentally measured Az; for the 



mono- and bi- layers and are indeed appreciably smaller 
than their geometrical heights. In particular, we find that 
an adsorbed multilayer corresponds to a vacuum layer 
with an height being about half the geometrical height of 
the multilayer. The increasing difference in geometrical 
height and apparent heights with the number of layers 
makes it increasingly difficult to image the layer. The 
larger Azi than Az; for Z > 1 is consistent with the be- 
havior of the calculated partial s, p, and d DOS around 
the CI atoms, as shown in Fig. [7| The interface layer 
has the largest conductivity because of the formation of 
a weak covalent bond between a CI atom in the inter- 
face layer and the Cu substrate results in a non- vanishing 
density of states in the gap region in the interface layer. 
Whereas for the other layers, this contribution to the den- 
sity of states are vanishingly small resulting in a smaller 
conductance and smaller apparent heights than for the 
interface layer. Note that the variation of the the work 
functions with the number of adsorbed layers contribute 
to Az/ but cannot explain the variation in Az/ with I. 
For instance, the decrease of the work function from the 
monolayer to the bilayer decreases the difference between 
Azi and Az2 because wave functions for electrons at the 
Fermi level decreases more slowly with decreasing work 
function. 

Finally, how the tunneling through a NaCl overlayer 
is affected by the presence of a band gap in bulk NaCl 
is obtained from the measured ^ spectrum. The char- 
acteristic features of this spectrum such as the strong 
reduction of ^ for voltages in the range from -4.5 to 2 
V is reproduced rather well by the simulated ^ spec- 
trum. This drastic reduction suggests a depletion of the 
LDOS in the corresponding energy region, which can be 
attributed to the absence of propagating states in the 
overlayer in this energy region and the formation of a 
band gap. However, our discussion in Section IlIII of the 
simulated ^ spectrum in Figs. showed that this sug- 
gestion is oversimplified. The onset at around y = 2 V is 
an LDOS effect whereas the onset &tV — -4.5 V is caused 
by the voltage-dependent barrier and field-emission2&. 



V. CONCLUSIONS 

In conclusion, we have carried out a combined scanning 
tunneling microscopy and density functional study of the 



surface topography, the apparent heights, work functions 
and the differential conductance of mono-, bi-, and tri- 
layers of NaCl on a stepped Cu(311) surface. The exper- 
imental results are well reproduced by the calculations, 
which also provide physical insights and interpretation of 
the results. The adsorption-induced changes of electronic 
and geometric structure of the NaCl overlayer are found 
to be significant. The observed large adsorption- induced 
changes in the work function are found to be associated 
with substantial charge rearrangements upon adsorption. 
The CI atoms were shown to be imaged as protrusions. 
In the STM images, also the different positions of the 
CI atoms relative to the underlying Cu(311) lattice in 
the p(3x2)-I and p(3x2)-II structures of the NaCl mono- 
layer were resolved. The poor conductivity of the layers 
were revealed by the apparent heights of steps being only 
half of the geometrical heights. Furthermore, ^ spectra 
revealed a large reduction of the tunneling conductance 
in a wide voltage region, resembling a band gap. How- 
ever, the simulated spectrum showed that only the onset 
at positive sample voltages may be viewed as a conduc- 
tion band edge, whereas the onset at negative voltages is 
caused by the drastic effect of the electric field from the 
tip on the tunneling barrier. 
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